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1. Introduction 



The concepts of scale invariance and scaling behavior are now increasingly 
applied outside their traditional domains of application, the physical sciences. 
Their application to financial markets, initiated by Mandelbrot [1,2] in the 
1960s, has experienced a regain of interest in the recent years, partly due to 
the abundance of high-frequency data sets and availability of computers for 
analyzing their statistical properties. 

This lecture is intended as an introduction and a brief review of current 
research in a field which is becoming increasingly popular in the theoretical 
physics community. We will try to show how the concepts of scale invariance 
and scaling behavior may be usefully applied in the framework of a statistical 
approach to the study of financial data, pointing out at the same time the 
limits of such an approach. 
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2. Statistical description of market data 



Quantitative research mainly focuses on "liquid" financial markets i.e. orga- 
nized markets where transactions are frequent and the number of actors is large. 
Typical examples are foreign exchange markets, organized futures markets and 
stock index markets and the market for large stocks. Prices are recorded sev- 
eral times a minute in such markets, creating mines of data to exploit. Such 
markets are complex systems with many degrees of freedom [3] , where many 
internal and external factors interact at each instant in order to fix the transac- 
tion price of financial assets. Various factors such as public policy, interest rates 
and economic conditions doubtlessly influence market behavior. However, the 
precise nature of their influence is not well known and given the complex nature 
of the pricing mechanism, simple deterministic models are unable to reproduce 
the properties observed in financial time series. Furthermore, although the de- 
tails of the price fixing mechanism- market microstructure - may be different 
from one market to the other, what is striking is the universality of some simple 
statistical properties of price fluctuations, prompting a unified approach to the 
study of different types of markets. 

As in the case of other types of complex systems with universal character- 
istics, a stochastic approach has proved to be fruitful in this case. The main 
object of study in this framework is the probability density function (PDF) of 
the increments at a given time scale T i.e. the probability distribution Pt of 
x(t + T) — x(t) where x(t) is the price of the asset at time t. This approach was 
inaugurated by Louis Bachelier who first introduced the idea that stock market 
prices behave as a random walk [4], and who considered Brownian motion as 
a candidate for modeling price fluctuations. 

Bachelier's model, applied to the logarithm of the prices (to ensure positivity 
of the price!), became very popular in the 1950s [5] and it is one of the main 
ingredients of the famous Black-Scholcs option pricing formula [6] . This model 
implies that the increments of asset returns (or asset prices) are independent 
identically distributed (iid) Gaussian variables. Indeed if one considers each 
price change as a sum of many small and independent random contributions 
from various market factors, the Central Limit Theorem suggests the Gaussian 
as a natural candidate. 

Figure 1. gives a typical example of the empirical distribution of the in- 
crements of asset prices, in this case the U.S. Dollar/ Deutschemark exchange 
rate, sampled every 5 minutes. The contrast with the Gaussian is striking: this 
'heavy tailed' or 'leptokurtic' character of the distribution of price changes, 
has been repeatedely observed in various market data and may be quantitavely 
measure by the kurtosis of the distribution Pt defined by : 



K = (A ^- A r )4 -3 (1) 

cr(T) 4 



where <j(T) is the variance of Atx = x(t+T)—x(t). The kurtosis is defined such 
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that k = for a Gaussian distribution, a positive value of k indicating a slow 
asymptotic decay of the PDF. The kurtosis of the increments of asset prices 
is far from its Gaussian value: typical values for T=5 minutes are: n ~ 74 
(US$/DM exchange rate futures), k ~ 60 (US$/Swiss Franc exchange rate 
futures), Hi — 16 (S&P500 index futures) [7, 8]. 

These observations imply that using a Gaussian PDF systematically under- 
estimates the probability of large price fluctuations, an issue of utmost impor- 
tance is financial risk management. 

Density of 5-minute increments 

USD/DEM exchange rate futures 



0.20 




30.0 



Fig. 1. — Probability density of 5 minute increments of USD/Deutschemark exchange 
rate futures. The lower curve is a gaussian with same mean and variance. 



3. Scale invariance and stable laws 

But why should stock price fluctuations have scale invariant properties in the first 
place? The answer lies in a generalized version of the Central Limit Theorem [9] : the 
distribution of the sum a large number of independent identically distributed random 
variables belongs to a family of distributions known as stable or stable Levy distribu- 
tions, characterized by their Fourier transform (characteristic function) [9, 10]: 

4>t,{z) =exp[-a M |^(l + i/3tan(/i7r/2)-^-)] (2) 

\ z \ 

The parameter /i, called the stability index, belongs to the interval ]0,2], /i = 2 
corresponding to the Gaussian distribution. /3 is a skewness parameter, (3 — for 
a symmetric distribution. Levy distributions are characterized by the property of 
being stable under convolution: in other words, the sum of two iid Levy-distributed 
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variables is also Levy distributed, with the same stability index li. More precisely, if 
Xi,i = l..n are n iid /i-stable random variables then the renormalized sum: 

has the same PDF as the JQ. In particular, the sum scales as N 1 ^ and not as yf~N 
which is the case in diffusive random walks. Levy distributions, in addition to their 
stability under convolution, have other interesting properties: except for the Gaussian 
(fj, = 2) all /U-stable PDFs have power-law tails with exponent 1 + fx: 

C 

M*) ~ -^T+JI ( 3 ) 
leading to an infinite variance and heavy tails. 

These observations led Mandelbrot to propose stable Levy distributions as candi- 
dates for the PDF of price changes of financial assets [1,2]. Mandelbrot observed 
that stable distributions offer a heavy tailed alternative to Bachelier's model while 
still enabling a justification of the model in terms of the Central Limit Theorem (see 
above). Furthermore, their stability under convolution gives rise to the scale invari- 
ance of the process: if appropriately rescaled, the increment at scale Nt will have the 
same distribution as the increment at scale r: 

The above relation means that the price process x(t) is self-similar with a self- 
similarity exponent which is the inverse of the stability index li. Self-similarity in 
market prices was first remarked by Mandelbrot in the 1960s in his seminal work 
cotton prices [1]. Recent studies have confirmed the presence of self-similarity and 
scale invariant properties in various markets: the Milan stock exchange [11], the 
S&P500 index [12], the CAC40 (Paris stock market) index [13] and foreign exchange 
markets [14, 15] and individual french stocks [16]. The value of li may be estimated 
for example by examining the scaling behavior of the probability of return to the 
origin Pt(0), which scales as j^jp, , which is indeed the method used in [12, 13]. The 
value of li found depends on the market considered but all values fall between 1.5 and 
1.7, n being lower for more volatile markets. 



4. Beyond scale invariance: truncated Levy nights 

Although the representation of price increments as Levy-stable iid variables accounts 
for their leptokurtic character as well as their self-similar properties, it has several 
drawbacks: careful comparison of the distributions at various time scales show that in 
fact it is not scale invariant [8, 15, 7]. The self-similarity properties described above 
do not hold at all time scales but only for short ones, typically less than a week. For 
example, if the price followed a Levy walk, different methods of estimating the self- 
similarity parameter would give the same value, namely 1/fj,. Yet in practice, data 
from different time resolutions give different values of the self-similarity exponent, 
between 0.5 (the value for a Brownian diffusion) and 0.6 (= 1/fi with n = 1.7). 
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Variance of price increments 



Scaling behavior of kurtosis 

S&P Index futures, 1991-95 
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Fig. 2. — Left: Scaling behavior of the variance of price increments for S&P 500 
index futures: the variance increases linearly with the time scale, a property which is 
consistent with the absence of significant linear correlation. Right: Scaling behavior 
of the kurtosis of price increments for S&P 500 index futures: the kurtosis decreases 
more slowly than in the case of iid increments where it decreases as 1/N. The solid 
line represents a fit with k(T) = k(V)/VT. 



Furthermore, the best fit by a stable Levy distribution [12] still overestimates the 
probability for extreme fluctuations: in other words, Levy distributions have tails that 
are 'too fat' compared to fluctuations of real prices. Tails of empirical distributions 
may be studied more closely by rank-ordering techniques [8, 17] (also called Zipf 
plots). They turn out to decay more slowly than a Gaussian but more quickly than 
any Levy-stable distribution: their decay is better described by an exponential tail. 
Last but not least, real price changes turn out to have finite variance [8, 15]. 

These apparently contradictory aspects- self-similarity at short time scales, break- 
down of scaling for longer time scales, truncated power law tails, finite variance- may 
be blended together into a consistent picture by a truncated Levy flight description 
[18, 8, 19]. Consider a random walk with increments distributed with an exponentially 
truncated Levy density p(x) [20]. p(x) may be characterized by its Fourier transform 
(characteristic function): 

= exp ( -a,[ (( a 2 + *Y /2 C0S 0* arctonff/a)) _ ^ 

COS (TTfl/Z) 

The above expression differs from the characteristic function of a stable law by 
the presence of a cut-off parameter a. 4>{z) is defined such that it has a regular 
Taylor expansion at the origin, assuring the finiteness of all moments. In terms of the 
variable x, p(x) decays exponentially for large values of x, which describes adequately 
the decay of the tails of the PDF of empirical price changes. The asymptotic behavior 
of (j> is given by 

00) ~ <?V0) 

i.e. for small values of x, p(x) behaves like a Levy-stable law of index /i. We can 
now proceed to describe the behavior of a random walk with iid increments having 
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a truncated Levy PDF p(x). For short time scales, p(x) behaves like a stable law of 
index /x: the short term dynamics of the process is therefore correctly described by 
a Levy flight which gives rise to the scale invariance with a self-similarity exponent 
H = [1, 2, 11, 12]. However, since the variance of the distribution is finite, the 

Gaussian version of the Central Limit Theorem applies and for large time scales 
T >• T* = kt, where k is the kurtosis of the distribution P T {x). Thus at large time 
resolutions the distribution Px will approach a Gaussian, which is consistent with 
the fact that the kurtosis tends to zero when T — > co (see figure). Between these 
two regimes -Levy flight for short times and Gaussian diffusion for long times- is a 
crossover regime which is characterized by the appearance of fluctuations of the order 
1/q. 




10* r 



10"' 



Sx 



Fig. 3. — Rank ordering of the increments of S&P 500 index futures. The cumula- 
tive distribution function of price increments is well represented by an exponentially 
truncated Levy distribution with finite variance. 



5. Correlation and dependence 

It is a well-known fact that price movements in liquid markets do not exhibit any 
significant autocorrelation: the autocorrelation function of the price changes 
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C(T] = Sx t 8x t+T - Sxt Sx t+T 
var(8xt) 



rapidly decays to zero in a few minutes: for T > 15 minutes it can be safely 
assumed to be zero for all practical purposes [19]. The absence of significant linear 
correlations in price increments and asset returns has been widely documented [21, 22] 
and often cited as support for the "Efficient Market Hypothesis" [23]. The absence 
of correlation is intuitively easy to understand: if price changes exhibit significant 
correlation, this correlation may be used to conceive a simple strategy with positive 
expected earnings; such strategies, termed arbitrage, will therefore tend to reduce 
correlations except for very short time scales, which represent the time the market 
takes to react to new information. This correlation time is typically several minutes 
for organized futures markets and even shorter for foreign exchange markets. 

The fast decay of the correlation function implies the additivity of variances: for 
uncorrelated variables, the variance of the sum is the sum of the variances. The 
absence of linear correlation is thus consistent with the linear increase of the variance 
with respect to time scale. 

However, the absence of serial correlation does not imply the independence of the 
increments [7] : for example the square or the absolute value of price changes exhibits 
slowly decaying serial correlations. Figure 5 displays the autocorrelation function 
g{T) of the square of the increments, defined as: 



_ ^ X t^ x t+T ^ x t ^ X t+T _ ^ X t^ x t+T ^ x t ^ X t+T 
var(Sx 2 ) /i4(r) — a(r) 4 



for S&P 500 Futures. The slow decay of g(T) (see figure) is well represented by a 
power law [7]: 



a = 0.37 ±0.037 5o = 0.08 

K 



This slow relaxation of the correlation function g indicated persistence in the scale 
of fluctuations. Other measures of the scale of fluctuations, such as the absolute 
value of the increments, exhibit the same type of persistence, a phenomenon which 
can be related to the "clustering of volatility", well known in the financial literature: 
a large price movement tends to be followed by another large price movement but not 
necessarily in the same direction. 
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Autocorrelation function of price changes 

USD/Yen exchange rate fit* data 1992-94 



Autocorrelation of square of price changes 

S&P 500 Index futures, 1991-95 



- Power law fit 
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Time lag : N = T/ 5 minutes 



Fig. 4. — Left: Autocorrelation function of price increments for the US Dollar/ Yen 
exchange rate. Right: The autocorrelation function g of square of price increments 
of S&P futures exhibits a slow decay, well represented by a power law with exponent 
a ~ 0.37. 

The persistent character of the scale of fluctuations may be quantitatively related to 
the anomalous scaling behavior i.e. slow decrease of the kurtosis as the time resolution 
increases. As a result, the convergence to the Gaussian is slowed down even more 
than in the case of a truncated Levy flight with iid increments. 

The presence of such subtle correlations points out the limits of a representation 
of market prices by a random walk model and the need to take into account finer 
effects due to nonlinear correlations and non-stationarity [24,7]. These effects are 
especially important in options markets, where the behavior of volatility -the scale of 
price fluctuations- is crucial in determining the value of the option. It has been shown 
[24] that the correlations present in the volatility -as represented by the function g 
defined above- are accurately reflected in market prices for options. 



6. Turbulence and finance 



Recently some attempts have been made to draw analogies between scaling properties 
of foreign exchange rates and turbulent flows [25], based on some similarities between 
their probability density functions. Indeed, the statistical approach to the problem 
is the same, price increments playing the same role as velocity increments in the case 
of turbulence and a formal analogy may be very tempting. However, in the light of 
the above remarks, it is clear that such an analogy cannot be pushed too far: there 
is an essential difference between the two phenomena, namely, the presence of strong 
correlations between velocity increments, leading to Kolmogorov's famous fc~ 5 / 3 law 
[26], which are absent in the case of price increments [19]. As remarked above, the 
autocorrelation function of price increments decreases rapidly to zero, resulting in a 
flat ("white-noise") spectrum (see figure below, taken from [19]). Translating into 
a power spectrum we obtain a dependence well approximated by oj~ 2 , which merely 
reflects the absence of correlation, as opposed to the lo~ 5 ^ 3 law for turbulent flows. 
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Fig. 5. — Fourier transform of the price autocorrelation function {x(t)x(t + r)) as 
a function of temporal frequency for the DEM-USD exchange rate (Oct 1991-Nov 
1994) (bottom curve) as compared to the A; -5 / 3 power- law spectrum observed for the 
spatial velocity fluctuations in turbulent flows (the data were recorded by Y. Gagne 
in a wind tunel experiment at R\ — 3050) (top curve). Inset: Fourier transform of 
the price change autocorrelation function (Ax(t)Ax(t + r)), which is completely flat 
(and would behave as the 1/3 —th power of frequency in turbulent flows). Units are 
arbitrary. 



7. Conclusion 



The above example illustrates that although methods of statistical physics and in 
particular the concepts of scaling and scale invariance may be successfully applied in 
the study of market price fluctuations, sometimes revealing unsuspected properties of 
financial time series, analogies with problems in physics should be handled with care 
in order to avoid erroneous conclusions. Keeping this point in mind, the scientific 
study of financial markets proves to be a fascinating subject in itself, in particular 
for the physicist whose theoretical tools may prove to be useful in uncovering new 
properties and mechanisms in financial data. 
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